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Abstract: Dark matter and dark energy are dominating components of the Universe. 
Their presence affects the course and results of processes, which are driven by the gravita¬ 
tional interaction. The objective of the paper was to examine the inffuence of the dark sector 
on the gravitational collapse of an electrically charged scalar field. A phantom scalar field 
was used as a model of dark energy in the system. Dark matter was modeled by a complex 
scalar field with a quartic potential, charged under a ?7(l)-gauge field. The dark com¬ 
ponents were coupled to the electrically charged scalar field via the exponential coupling 
and the gauge field-Maxwell field kinetic mixing, respectively. Complete non-linear sim¬ 
ulations of the investigated process were performed. They were conducted from regular 
initial data to the end state, which was the matter dispersal or a singularity formation 
in a spacetime. During the collapse in the presence of dark energy dynamical wormholes 
and naked singularities were formed in emerging spacetimes. The wormhole throats were 
stabilized by the violation of the null energy condition, which occurred due to a significant 
increase of a value of the phantom scalar field function in its vicinity. The square of mass 
parameter of the dark matter scalar field potential controlled the formation of a Cauchy 
horizon or wormhole throats in the spacetime. The joint impact of dark energy and dark 
matter on the examined process indicated that the former decides what type of an object 
forms, while the latter controls the amount of time needed for the object to form. Addi¬ 
tionally, the dark sector suppresses the natural tendency of an electrically charged scalar 
field to form a dynamical Reissner-Nordstrom spacetime during the gravitational collapse. 
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1 Introduction 

Black holes which exist in the Universe are rotating and electrically neutral, due to the pres¬ 
ence of an interstellar medium. They are formed during gravitational collapse of a star, 
whose mass exceeds the Tolman-Oppenheimer-Volkoff limit [1, 2]. One of theoretical 
toy models of this process, which allows to describe the resulting spacetime structures 
and properties of the emerging objects, is a dynamical evolution of an electrically charged 
scalar field [3, 4]. A considerable advantage of its investigations is that they allow to describe 
physics of both outer and inner regions of the arising objects. 

Dynamical gravitational collapse attracts scientific interest also because the structures 
of dynamically formed spacetimes differ significantly from their non-dynamical counter¬ 
parts. This means that the knowledge of static solutions, which describe a particular 
system, is insufficient to characterize the results of dynamical processes, which take place 
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within it. For this reason, the evolution of a massless electrically charged scalar field under 
the influence of gravity was examined [5, 6]. It turned out that the dynamical Reissner- 
Nordstrom spacetime is substantially different from the static one due to the existence 
of a central spacelike singularity and a null singularity along the Cauchy horizon, which 
results from the mass inflation phenomenon [7]. The effect of pair creation during the col¬ 
lapse and the subsequent evaporation of a black hole were also considered [3, 4, 8, 9]. 
The collapse of a Brans-Dicke field was elaborated in [10-13], while the results of the gravi¬ 
tational evolution in Einstein-Maxwell-dilaton theory with standard and phantom couplings 
of the involved fields were described in [14-16]. 

According to current cosmic microwave background (CMB) measurements performed 
by Planck, the amount of dark matter and dark energy in the Universe is about 27% 
and 69%, respectively [17]. Such high abundance of the dark components, both of which 
interact gravitationally, undoubtedly affects the course and results of the gravitational col¬ 
lapse, which was one of the motivations for the analyses presented in the paper. The con¬ 
ducted research concerning the dark sector influence on the relatively small-scale process 
supplement vast investigations on its significance during the cosmological evolution. 

Until now, the collapse in the presence of the dark sector was studied mainly in the con¬ 
text of the large scale structure formation and the course of a few astrophysical processes. 
The evolution of dark matter concentration in the form of a spherically symmetric dust 
cloud of a finite radius was examined on a background spacetime containing dark energy 
modeled by isotropic and anisotropic fluids, as well as the Chaplygin gas [18-20]. The struc¬ 
ture formation during the evolution of quintessence dark energy induced by the collapse 
of dark matter halos was described in [21]. The spherical collapse model of the structure for¬ 
mation on cosmological scales in the presence of the dark sector was also considered [22, 23]. 
The evolution was followed on a background of a flat, homogeneous and isotropic universe, 
while the matter content was modeled by a multicomponent fluid. The clustering process 
was investigated during the non-linear spherical ’top hat’ collapse of the generalized Chap¬ 
lygin gas, which unifies dark energy and dark matter into a single matter component [24]. 

Recently, the role of the dark sector in gravitational collapse has been studied from 
the astrophysical perspective. The issue of a neutron star collapse due to capture and sed¬ 
imentation of dark matter within its core was considered as a solution to two interesting 
observational puzzles [25]. One of them are fast radio bursts, which are probably extragalac- 
tic, transient bright radio pulses of an unexplained origin [26]. The second is a missing pulsar 
problem, which is related to the non-detection of pulsars within the inner core of the Galac¬ 
tic Center of the radius close to 10 parsecs [27]. The early stage and the long-term phase 
of the formation of supermassive black holes progenitors via a direct gas collapse in dark 
matter halos at high redshifts was studied in [28, 29]. 

The above-mentioned analyses were mostly focused on the behavior of matter during its 
gravitational collapse. The research was performed within the framework of fluid dynamics 
on either fixed or evolving background spacetimes. The current studies, by contrast, involve 
examining dynamical matter-geometry systems with the matter content of spacetime de¬ 
scribed in terms of classical field theory. To our knowledge, such a methodology has not yet 
been used in the analyses of the role of the dark sector in gravitational collapse. It allowed 
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US not only to investigate the matter behavior, but also to follow the geometry dynamics 
observed during the evolution. The latter was reflected in the obtained spacetime struc¬ 
tures. Due to the fact that they resulted from the fully non-linear computations conducted 
within the whole spacetime region subjected to the collapse dynamics, they were not lim¬ 
ited by the capabilities of the previously employed local techniques [18, 19] and allowed 
presenting global characteristics of the forming spacetimes. Their interpretation enabled 
us to describe causal relations between distinct spacetime areas and to test the geometrical 
properties of both outer and inner regions of the emerging objects. 

The first indirect evidence suggesting the existence of an electromagnetically non-inte¬ 
racting unknown matter component of the Universe came from observations of the galactic 
rotation curves [30, 31]. Since then, the concept of dark matter has put down its roots 
in modern cosmology. Apart from the aforementioned experiments, the gravitational in¬ 
teraction of dark matter with the remaining content of the Universe may be inferred from 
the observations of CMB, baryon acoustic oscillations (BAO), the rate of the large scale 
structures formation [17, 32] or galactic collisions, like the Bullet Cluster [33]. These obser¬ 
vations also provide information of the dark matter abundance and distribution in the Uni¬ 
verse, as well as its non-gravitational nature, i.e., possible interaction types and the lifetime 
of the dark matter particle. Unfortunately, they are insufficient to fully identify the particle 
nature of dark matter, which is a dominant matter component of the Universe. However, 
the non-gravitational interactions of dark matter disfavoring some of the proposed exten¬ 
sions of the Standard Model were revealed recently [34, 35]. 

There is yet another branch of modern physics that is immensely interested in the topic, 
namely particle physics. The Standard Model of particle physics, despite being very suc¬ 
cessful in describing the visible matter, lacks any dark matter candidates. From this per¬ 
spective the quest of gathering and analyzing information about dark matter is the quest 
of understanding the beyond-Standard Model physics. This situation and an enormous 
technological advancement triggered in the last few decades made the subject of dark mat¬ 
ter one of the most interesting in modern physics. Astrophysical observations (other than 
these mentioned above) provide hints that dark matter may be responsible for the y-ray 
bursts from galactic centers [36-38] and keV photons from galactic clusters [39, 40]. The ter¬ 
restrial experiments in supercolliders such as the Large Hadron Collider (LHC) also strive 
to identify dark matter or at least constrain its models [41], for a current list of the dark 
matter models proposed to search for in the LHC run-2 data see [42]. Apart from the above 
two sources of the information on dark matter, a new one appeared recently. It is con¬ 
cerned with the theoretical investigation of the role of dark matter in the missing pulsar 
problem [25, 43] and its impact on the composition and lifetime of the early generation 
of stars [44]. 

As a dark matter model we used one of the viable, from the particle physics per¬ 
spective, models of the dark matter sector. It is composed of a complex scalar field with 
a quartic self-interaction, which is charged under an additional Abelian gauge field. Its cou¬ 
pling to the standard matter sector (represented by an electrically charged scalar field) 
was implemented via a kinetic mixing between the introduced gauge field and the elec¬ 
tromagnetic field [45]. A new experiment constraining such a coupling is currently under 
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preparation [46], while some previous results can be found, e.g., in [47]. The considered 
model describes one or two dark matter candidates, depending on the vacuum expectation 
value (vev) of the complex scalar field. If the vev is non-zero, the dark matter candi¬ 
date is the massive gauge boson mentioned earlier (in the particle physics literature often 
called Z' or dark photon). When the scalar does not possess the vev the candidate can 
be either the scalar or the gauge boson. A slightly modified version of this model was 
analyzed in the context of dark matter phenomenology [48, 49] and the GeV-scale y-ray 
excess in the Galactic Genter [50]. On the other hand, a new theoretical tool for detecting 
the influence of the dark matter sector on characteristic quantities of holographic supercon¬ 
ductors was proposed [51-53]. It happens that the coupling constant of the dark matter 
sector to the ordinary Maxwell field imprints its existence during holographic phase tran¬ 
sitions and these phenomena may constitute a possible way of understanding the nature 
of the dark matter sector. 

Dark energy is a notion assigned within the ACDM cosmological model to the content 
of the Universe, which remains after determining the amount of matter detectable either 
directly (photons, neutrinos and baryonic matter) or indirectly (dark matter). It is re¬ 
sponsible for the present-day accelerated expansion of the Universe [54-57]. Experimental 
restrictions on dark energy models are imposed by observations of supernovae la, combined 
with the GMB and BAG measurements [58]. The obtained constraints allow of the existence 
of dark energy in the form of phantom matter [59]. Its pressure P and density p satisfy 
the relation P < —p^ which means that the value of its barotropic index w = Pp~^ is less 
than —1. Matter of this type is coupled repulsively to gravity and thus exerts negative 
pressure and causes cosmic acceleration. Such a phantom scalar field coupled exponentially 
to the electrically charged scalar field was chosen as a dark energy model in our research. 

Investigating a dynamical evolution in the presence of phantom fields is also justified 
due to a wide variety of non-dynamical solutions to Einstein equations involving such fields. 
New classes of spherically symmetric solutions of Einstein-Maxwell-dilaton theory with 
a phantom coupling were considered in [60, 61], while static multicentered solutions in this 
theory were discussed in [62]. A vast array of diversified and complex solutions opens new 
possibilities of finding their interesting dynamical counterparts. 

The concept of a static wormhole was introduced by Wheeler [63] and extended by Mor¬ 
ris and Thorne [64], who invented a class of potentially traversable wormholes. The worm- 
hole throat is kept open by the presence of exotic matter, whose energy-momentum tensor 
violates the energy conditions. A scalar field coupled to gravity in a phantom manner 
is a good candidate for such a type of matter, which supports wormhole traversability [65]. 
The quoted works initiated studies of wormholes as topological bridges connecting two dis¬ 
tinct regions of spacetime. There were proposed various models such as thin shells in cos¬ 
mological background [66], in Brans-Dicke gravity [67], in the low-energy limit of the string 
theory [68], as well as wormholes supported by scalar phantom fields [69, 70], tachyonic mat¬ 
ter [71] or non-linear electrodynamics [72]. Recently, to circumvent the problem of exotic 
matter, the modified-gravity theories were considered as suitable candidates to obtain sta¬ 
ble traversable wormholes [73-76]. It is worth emphasizing that these researches deal with 
constructing eternal, i.e., non-dynamical wormholes. There also exists a research branch. 
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which is concerned with physics of evolving wormholes [77-79]. The aim of our studies 
was to test the dynamical formation of such objects during fully non-linear investigations 
of the gravitational collapse. 

Because of the ambiguous nature of wormholes in the view of astrophysical observations, 
a considerable resurgence of interests in their physics was observed recently. Experimental 
studies of phenomena and processes which result from the presence of a black hole in space- 
time are based on the effects related to its external gravitational field. Wormholes possess 
the same gravitational properties as black holes, which are characteristic for them from 
the observational viewpoint [80]. For this reason, a wormhole is an alternative to a black 
hole and hence various possibilities of the experimental distinction between these two types 
of objects are being analyzed [81, 82]. 

Apart from the efforts to observationally distinguish wormholes from black holes, the is¬ 
sue of experimental distinction between naked singularities and black holes has also been 
pursued. These two types of objects could be differentiated through their gravitational 
lensing features [83] , such as a number of generated images of a particular light source [84] , 
their orientation in space, total magnifications and time delays [85, 86]. Another possi¬ 
bility is investigating the properties of accretion disks, which form around naked singu¬ 
larities and black holes, as their luminosities and angular velocities of the particles within 
them depend on the type of the central object [87, 88]. The above findings indicated that 
the observational consequences of the existence of naked singularities and black holes differ. 
It opens a possibility of testing the cosmic censorship hypothesis experimentally through 
astronomical measurements [89]. 

Astrophysical objects present in the Universe are dynamical, as they constantly undergo 
changes due to interactions with the surrounding environment. Thus one of the elements 
necessary to achieve the aforementioned goal, that is observational determination of a true 
nature of the astrophysical objects with strong gravitational fields, is to investigate the for¬ 
mation and behavior of black holes, wormholes and naked singularities during dynamical 
evolutions. 

The leading research objective of the current studies was to test the inffuence of dark 
energy and dark matter on the collapse of an electrically charged scalar field. The role 
of the dark sector during the analyzed process was investigated through the prism of the form¬ 
ing spacetime structures, the behavior of fields within them and properties of the emerging 
dynamical objects. The paper is organized as follows. In section 2 we present the theoreti¬ 
cal model, which enables us to study the dynamical gravitational collapse of an electrically 
charged scalar field accompanied by the dark sector. Section 3 contains basic informa¬ 
tion on solving the derived equations of motion and particulars of the results presentation. 
In sections 4 and 5 we discuss the course and outcomes of the investigated process run¬ 
ning in the presence of dark energy and dark matter, respectively. Section 6 is devoted 
to studies of the joint inffuence of these two components on the evolution. The summary 
of the obtained results and prospects of further research are placed in section 7. Comments 
on numerical methods used for solving the equations of motion and the code accuracy checks 
can be found in appendix A. 
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2 Theoretical model of the evolution 


The model used to investigate the dynamical collapse of interest consists of three parts, 
whose construction is based on self-interacting scalar fields. The first one refers to a complex 
scalar field ^ coupled with the Maxwell field A^. The second one is a phantom scalar field 0, 
which represents dark energy in the system. Its exponential coupling to the previous field 
is consistent with the low-energy string theory regime. The third component is a massive 
complex scalar field x with a quartic potential coupled with a C/(l)-gauge field which 
is coupled to the Maxwell field. It serves as dark matter in the system. Since the model 
incorporates notions taken from string theory, it is convenient to write the most general 
form of the action in the string frame 




/ 




-2ct> 


R 


-2^(^V(/)^ +e^^^tsF d Cdm^ , (2.1) 


where the constant ^ determines the nature of the coupling between 0 and gravity. It is set 
as equal to —1, which means that the field is phantom [61, 62]. The constant a char¬ 
acterizes the coupling between the scalar field (j) and the electrically charged one. Since 
the low-energy limit of the string theory is considered, we set a as equal to —1 [90]. This 
value was incorporated in the equations below. The dark energy and dark matter sectors 
are chosen to be independent, so dark matter represented by the Lagrangian jCdm is not 
coupled to the phantom scalar field via the exponential coupling The presented 

model enables us to study the impact of dark components on the evolution both separately 
and jointly. Throughout the computations, we used the geometrized units system, in which 
SttG = c = 1. 

The Lagrangian of the electrically charged scalar field d’ is given by the expression 

CsF = -\d^^ (d^^)* - (2.2) 


where F^(j is the strength tensor of the Maxwell field. The covariant derivative has the form 
Dp = Vp + ieAp^ where e is the electric coupling constant, Ap is the four-potential and i de¬ 
notes the imaginary unit. 

The dark matter Lagrangian is the following: 

Cdm = -D^x {FxJ - - V (IxP) • (2.3) 


The covariant derivative is Dp = Vp + iePp^ while Pp is the four-potential of a U{1)- 
gauge field, e is its coupling constant with x and Bpcr = dp Pa — daPp- The scalar field 
potential is given by "L (|xP) = ^|xl^ + The constants m? and \dm are a square 

of the mass parameter of the scalar field and its quartic self-interaction coupling constant, 
respectively. The square of the mass parameter can be either positive or negative, when it is 
negative the scalar field possesses a non-zero vacuum expectation value. On a side note, 
due to the system of units, which was used in our calculations, it is impossible to compare 
the parameter m? with the particles masses directly. The perturbativity of the quantized 
theory requires that \dm is in general smaller than dvr. In the considered theory it is 
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restricted to values not exceeding 0.2 [48]. The parameter ajjM stands for the kinetic 
mixing coupling constant, which controls the strength of the mixing between the photon 
and the dark photon. Taking into account various experimental scenarios, its upper value 
is constrained to be not larger than 10“^ [91]. 

The searches for the dark photon have been conducted with the use of diverse ex¬ 
perimental techniques. The current constraints on its mass and its coupling constant with 
photon come from a collection of experiments, which are adapted to inspecting their various 
value ranges. The respective upper limits obtained during the dish antenna experiments 
are 3.1 eV and 6-10“^^ [92]. Studying the supernovae excess cooling covers the a^M values 
from 10“^^ to 10“^ and masses within the range 10“^ — 1 GeV [93]. The beam dump exper¬ 
iments [94, 95], searching for lepton jets in pp collisions [96], investigating heavy neutrino 
decays [97] and analyses of the distortion of the CMB [98] allowed examining the kinetic 
mixing parameter from within the range 10“^ — 10“^ for the dark photon mass ranges 
0.03 — 0.63 GeV, 0.4 — 1.1 GeV, 1 — 500 MeV and 10“^^ — 10“^ eV, respectively. Optical 
experiments dealt with the dark photon masses between 10“^ and 10“^ eV and the values 
of the mixing parameter greater than 10“^ [99]. The upper limit of a^M was established 
as equal to 1.7 • 10“^ for the mass range 5 — 470 MeV in experiments focused on ana¬ 
lyzing meson decays [100-102]. Investigating e^e“ collisions allowed searching the mass 
range 0.02 — 10.2 GeV [47], while analyses of the e~ and p scattering on nuclei focused 
on 175 — 550 MeV range [103-105], both with the values of the coupling constant with 
photon of the order of 10“^ — 10“^. 

The variation of the action (2.1) with respect to adequate fields, namely phantom 
scalar 0, Maxwell complex scalar fields and x? as well as the gauge field leads 
to the following set of equations of motion: 


= 0, (2.4) 

i - ie^ = 0, (2.5) 

+ ieA^ ( 2 V /3 + ieAf^) ij; + ieW 13 A^'ijj = 0, (2.6) 
VX* - ieA>^ (2V/3 - ieAf^) V’* - ieV/^A^^P* = 0, (2.7) 

X - e^PpP^x + ixP^) + Pp (V^x) - - ^^IxlX = 0, (2.8) 

1 


VX* - e‘^PpP^X* - 


m * Xdm I 
^X — 


IxlX* = 0 , (2.9) 


- 2e^Pnxr - *e(xWx* - X*Wx) + = 0.(2.10) 


During the derivation of the above evolution equations the string frame was converted into 
the Einstein frame in order to get manageable equations of motion. The metrics in these 
two frames are related via the conformal transformation 


g^lu = e (2-11) 

where and g^y denote metric tensors in the Einstein and string frames, respectively [106]. 
The transformation between frames (2.11) preserves causality. Since the interpretation 
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of the obtained results is based mainly on investigating notions related to causal structures 
of emerging spacetimes, we do not expect the transformation to influence the ultimate 
conclusions, which can thus be regarded as physically relevant. From now on, all quantities 
are written in the Einstein frame. 

The issue of choosing an adequate conformal frame to describe physics of a particu¬ 
lar system has been so far addressed on numerous occasions in research related primarily 
to cosmology and particle physics (see, e.g., [107-110]). A general conclusion on the equiv¬ 
alence of conformal frames at the classical level is that although the frames are physically 
equivalent, the interpretations of results obtained within them may differ [111]. However, 
as was stated above, this drawback does not refer to our case, as the results analysis is based 
mainly on describing causal properties of spacetimes. 

The most convincing arguments on physical relevance of a particular frame could be ob¬ 
tained by comparing theoretical and experimental data. Some attempts in this respect have 
been made recently in relation to cosmological observations [111-113]. However, they are 
not yet complete enough to resolve the issue of the conformal frames equivalence satisfac¬ 
torily and hence the problem of its experimental confirmation still remains open. 

The Einstein equations derived by varying the action (2.1) with respect to gravitational 
field complement the above set of equations, which describes the examined dynamical sys¬ 
tem. For studying the evolutions of interest, the double null spherically symmetric line 
element [114] is selected 

ds‘^ = —a{u^ v^dudv -\- v)dVt^^ ( 2 . 12 ) 

where u and v are retarded and advanced time null coordinates, respectively, and = 
dQ^ -h sin^ Qd^^ is the line element of the unit sphere, where 0 and $ are angular coordi¬ 
nates. Such a coordinate choice determines the spacetime foliation for conducting computa¬ 
tions, which is 2+2 [115]. The double null formalism was successfully employed in dynamical 
gravitational collapse investigations, e.g., [3-6, 8, 9], because it enables to follow the evolu¬ 
tion from approximately past null infinity, through the formation of horizons up to the final 
central singularity in the case of singular spacetimes. 

Regarding spherical symmetry, the only non-vanishing components of the field tensors 
are Fy^, Buv and Byy. Due to the gauge freedom Ay AyF'Vu^' and Py PuA'VuO", 
where 0 ' = J Aydv and 9" — J Pydv^ the only non-zero four-vector components are Ay 
and Pu- They are functions of retarded and advanced time. 

The equation of motion for the phantom scalar held (2.4) in the chosen coordinate 
system is provided by 

1 O^o? 

r,u(t>,v + + r(t)^uv + 

where we set 

4 

Q = 2^^. (2.14) 

Q is a function of retarded and advanced time, which corresponds to electric charge within 
a sphere of a radius r{u^v)^ on a spacelike hypersurface containing the point {u^v). Partial 




derivatives with respect to the null coordinates are marked as and Concerning the as¬ 
sumed line element (2.12) and the definition of electric charge (2.14), the 'u-component 
of Maxwell equations (2.5) can be separated into two first-order differential equations. 
The first one governs the evolution of the only non-zero component of the four-vector 
of the Maxwell field 

Oa^ 

1:^ = 0’ (2.15) 

while the second one describes the dynamical behavior of Q. Namely, one has 

- V’V’t;) + ^oiDMe^'^Tv = 0, (2.16) 

4 ’ O 

where we denoted by T the expression 

T = —Pu,v, (2.17) 

being related to the charge associated with the field within a sphere of a radius r{u^v) 

on a specific spacelike hypersurface containing The relations for the complex scalar 

field (2.6)-(2.7) imply 




uy iCT Aylp^y “h iCT^y Ayljj “h 
- ierAu^'^y - ier^yAy^"' - 


ieQa 

4r 

ieQo? 

4r 


2 

= 0 , 

■7/^* = 0. 


(2.18) 

(2.19) 


The equations of the scalar field y, its complex conjugate x* and the related ?7(l)-gauge 
field are given by 


+ r,uX,v + r,vX,u + ^iexrPu,v + ^ierPuX,v + ^P„x + 

I 1 2 2 I ^DM 2 I |2 n 

+ -m a rx+ “^Ixl X = 0) 
rx*uv + r,uX% + r,vX% - \iex*rPu,v - herPuX% - her^^PuX* + 


I ^ 2 2 ^1 

+-m a rx + 


2 

X£)M 2 


a ^'Ixl X* = 0, 

o o 

2ier^ {xx*v - X*X,v) + iz^dmQ^v = 0, 

Ta^ 

p _— 0 

2r2 ~ 


The stress-energy tensor for the considered theory is the following: 

^dtiuDp^p ^ D^'tp + 

F^X ( t ^ i ^ x ) + ( t ^ mx ) Di'X + BfipFif + oiDMB^pF^^ + 


4' 

1 

+ 2 


-dlJ-u 


Dpx (d^xJ + \bp„B^- + + V (1x0 


( 2 . 20 ) 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


(2.24) 
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Its non-vanishing components calculated in double null coordinates are of the forms 


Tuu — 2^0,n + 2 


ip,u'tp*u + ieAu (V’V’.u - V’*V’,'u) + 
+X,uX% + iePu {xx*u - X*X,u) + e^PuXX*, 


+ 


+ x,vx%, 


_ 

-‘-uv — ^ 


+ 




2r4 4r4 


+ O^dM 


4^4 

1 


A i 


TQa^ a 


Tee = + e-^^% + - ^ 

TO 2r^ 


4 V 2 

+ ieAu {'Ip'lp^ - '4’*'4’,v) 


+ 


TQ 


X,uX*v + X,vX*u + i^Pu {xx*v - x*x,v) 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


Combining the adequate components of the Einstein tensor resulting from the metric (2.12) 
and the above stress-energy tensor components, the Einstein equations of the gravitational 
field are obtained 

r 




2 a .,r . 


a 


^ 


+ 9 X,uX*u + 'i^Pu {xx*u ~ X*X,u) + PuXX* 

2 


T T 

— 'r.7;7; = + l^X.vX^V') 


(2.29) 

(2.30) 


a , r^uTv , 

^-h r^uv — e 

4r r 


_ ^-20^ 


+ 


4r3 


+ 


8r3 


+ C^DM 


TQo? ro? (m? o 

^ + - ^lxf + 


^DM 


ixr 


^ _Lp}L = (, ^ i(t),u4>,v + ^ + ieAu i'tp'ipl - 'tp*'ip,v) 

4 


4r^ 


(2.31) 

+ 


8r4 


TQa‘^ 

+ OIDM —:r^i- 1 - 


8r4 


1 

+ 2 


X,uX*v + X*uX,v + iePu ixxl - x*x,v) 


(2.32) 


They complement the preceding relations in order to obtain the complete set of equations 
of motion for the examined system, which are (2.13)-(2.32), excluding the relations (2.14), 
(2.17) and (2.24), which determine the physical quantities Q, T and 

In order to make the solution of the obtained equations of motion attainable, one in¬ 
troduces a set of auxiliary variables 


(2.33) 


c 

®,u 

d — 


/ 

= ^5 

II 


a 


a 


k 

= 4>, 

X — 

0,115 

y 

II 


s 

= 'tp, 

p = 

0,115 

q 

II 

3 

II 

h 

= X, 

w = 

X,n5 

z 

= X,v, 

II 


and the quantities denoted by 

X - “T f 

A = ^ + /5, 


l^ = fq + 9P, K = gw + fz, 


(2.34) 
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which make it possible to rewrite the second-order differential equations (2.13), (2.18)- 
(2.21) and (2.29)-(2.32) as first-order ones. Moreover, the real fields t/^i, '02, Xi X2 
are introduced instead of conjugate fields 0, 0*, x X* according to 0 = 0i -h i02, 
0* = 01 — '^02, X = Xi + ^X 2 and X* = Xi ~ '^X 2 - These relations result in 


S ^ Si+ iS 2 , 

h = hi + i/i2, 
/i = /ii -h i/i2, 
K — Ki + m2, 


p = pi+ ip2, 

W = Wi + iW 2 ^ 

Ml = /gi + 

= fzi + gwi, 


q = qi+iq 2 , 

Z = Zl+ iZ 2 , 

M 2 = fq 2 +gP 2 , 
1^2 = fZ 2 + gW 2 . 


(2.35) 


The final system of equations of motion, which governs the investigated evolution yields 


PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

PI 


5, 


a,u = ac, 

Oj — ad^ 

r,u = /, 
r,v = g, 

Sl(2),u = Mi(2)) 

■ 81 ( 2 ),= ^1(2)) 

^1(2),u = ll’l(2)) 
hl{ 2 ),v = ^1(2)) 

/,„ = 2 cf -rx"^ Pi + P 2 + 2e/3 {sip 2 - S 2 M 1 ) + e^/3^ (sf + s^) 


+ 


1 

’ 2 ^ 


w 


I + + 2e7 (Miru2 - h2rui) + (hf + h^) 


E2 : = 2dg - ry‘^ - ^ {qi + q^) - ^ {zi + z^) , 


A 


- 2 kQ‘^a^ 


E2> : g,u — f,v — 1“ e „ 

rp ZlpO 


+ 


T^a? 


TQa^ 

+ ^ Q—h 


8r^ 


+ - 


ra 


Y ('‘l + '•M + ^ ('•1 + hE 


P4 : d „ = c ^ = A _ ajy _ 1 + p^q^ + (sig2 - S 2 Q 1 ) 


— e 


- 2 fce« 


2^2 


2r^ 


+ 


8 r^ 


o-DM- 


TQa^ 

8^ 


WlZl - W2Z2 + 67 (^ 22:1 - hlZ2) , 

Qo? 


(Re) ■ 


qi,u = Pl,v = - — + 6/3^2 + 652/3- + eS 2 , o , 

' ' p p llp^ 


\lm) 


M2 „ o9 Qa^ 

q2,u =P2,v = - epqi - 651/3 - 651 —^, 

^ ^ /y* ry* /yry^Zj 

Ml: /3. = 


2r2 ’ 


M2: Q,^ = 

OiDM ~2„2k 


1 


2ye + 2®^ ® (®192 - S2gi) + 


+ 


{zih 2 - Z 2 hi) 


2 \ “1 
^ 2 k ^DM 

^ 16 


(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


(2.46) 

(2.47) 

(2.48) 

(2.49) 

(2.50) 

(2.51) 
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H, 


(Re) • 




K.I ^ Ta^ 1^ 1^ g 

zi^u = =-h eh 2 —^ + -62:27 + 06/127- + 

r 8r^ 2 2 r 

1 2 7 

--a hi 


rn^ + Ai:)M (/?'i + hf) 


(Im) 


K2 Ta^ 1 , I g 

Z2,u = = - ehi—y - -ezi-f - -ehij- + 

r 8r^ 2 2 r 

— ^a^/i2 


+ \]jM (j^i ^ 2 ) 


Cl : = 


Ta^ 


C2 : T^y = 2er^ {hiZ2 — /i2^i) — 


Q.DM 


OiDM ~ 2 2 k ( u u \ 

H—^er e {zih 2 - ^ 2 ^ 1 ) 


1 — e 


2 klDM_ 
16 


2y<3 + (siQ2 - S2yi) + 

-1 


(2.52) 


(2.53) 

(2.54) 


(2.55) 


3 Details of computer simulations and results analysis 

Because of its complexity, the system of the obtained differential equations (2.36)-(2.55) 
needs to be solved numerically. The details of the numerical code and performed tests are 
presented in appendix A. 

The evolution equations were solved in the region of the (u?i)-plane, which is shown 
on the background of a dynamical Reissner-Nordstrom spacetime [6] in figure 1. In all 
conducted simulations it was confined to 0 ^ u ^ 7.5 and 0 ^ ^ 7.5. The only arbitrary 

input data of the computations were profiles of the evolving fields, posed on the initial null 
hypersurface denoted as ?/ = 0. The initial profile of the phantom scalar field was Gaussian 

C2 ; , 

while the complex fields were modeled by the trigonometric profile of the following form: 

V’(or x) = Psior Ph) ■ sin^ ( tt— 

V Vf 

The above profiles were selected so that they described the behavior of real and complex 
scalar fields properly [5, 6, 116, 117]. The profiles were treated as one-parameter families 
with amplitudes Ps and ph as free family parameters. The amplitudes are indicators 
of the strength of the gravitational self-interaction of the particular field [118]. The existence 
of a specific matter type in the examined system was guaranteed by the non-zero value 
of the respective amplitude. Three systems were investigated, i.e., 

• electrically charged scalar field - dark energy (SF-DE) with ps ^ 0, Pk ^ 0 
and Ph = 0, 

• electrically charged scalar field - dark matter (SF-DM) with ps ^ 0^ Pk = 0 
and Ph ^ 0, 


, 2 v\ . f 2 v 

cos I TT- -h I cos TT-h 0 

VfJ \ Vf 


(3.2) 



(3.1) 
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Figure 1. The computational domain (marked gray) 
on the background of the Carter-Penrose diagram of the dy¬ 
namical Reissner-Nordstrom spacetime. The central singular¬ 
ity along r = 0, the event and Cauchy horizons are denoted 
as iS, £% and CH, respectively. and are null and time¬ 
like infinities, while is a spacelike infinity. 


• electrically charged scalar field - dark energy - dark matter (SF-DE-DM) with 
Ps 7^0, Pk 7^0 and ph 7 ^ 0. 

The remaining constants were arbitrarily chosen and invariable during computations. 
To be precise, we set ci = 1.3, C 2 = 0.21 and the parameter determining the amount 
of initial charge 5 — The maximum value of advanced time was Vf = 7.5. The choice 
of initial conditions is representative for the conducted evolutions, because their outcomes 
are independent of the types of profiles provided that they are regular, i.e., they result 
in a regular spacetime slice at the initial null hypersurface. This condition is fulfilled 
by the selected profiles (3.1) and (3.2). 

Provided that the value of the electric coupling constant is non-zero, it does not affect 
the results of the collapse [14], which was confirmed for the investigated cases. For this 
reason, it was set as equal to e = 0.5 in all evolutions. A similar property was observed 
for the coupling constant e, so it was also kept constant and equal to 0.5 during simulations. 
The cases of vanishing e and e were not examined, because they refer to uncharged complex 
scalar fields, whose behavior during the collapse is beyond the scope of the current research. 

The dynamical spacetime structures resulting from the considered evolutions will be 
presented on Penrose diagrams. These diagrams contain contours of r = const, lines plotted 
in the ('L’i/)-plane. The outermost thick line refers to r = 0, which is non-singular when 
coinciding with the u — v line and singular in the remaining part. The lines indicating 
the vanishing expansion 


9- = -r ■ 


(3.3) 


r 


in the respective null directions i — u.^v will be also presented on the diagrams and their 
role in a spacetime will be commented in each case separately. They will be denoted 
as — 0 and r^u — 0 and marked on the diagrams as red and blue solid lines, respectively. 
The presentation of Penrose diagrams will be limited to spacetime regions, which are crucial 
for the conducted results analysis. 

The aim of the paper is to present the influence of the dark sector on the collapse 
of an electrically charged scalar field. In figure 2 we present the Penrose diagrams of space- 
times formed during the gravitational evolution of the field without any additional factors 
as a reference for further analyses. For the initial amplitude value ps — 0.2 the emerg¬ 
ing curved spacetime is non-singular, because the evolving field disperses towards infinity 
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(a) 


(b) 




Figure 2. (color online) Spacetime structures emerging from the collapse of an electrically charged 
scalar field of the amplitude (a) Ps = 0.2 and (b) Ps = 0.6. 


(figure 2(a)). The process with ps = 0.6 results in a dynamical Reissner-Nordstrom space- 
time [6], shown in figure 2(b). 

The obtained Reissner-Nordstrom spacetime consists of two parts. The first one covers 
the dynamical region of spacetime and spreads up to about v = 7.05. The second one refers 
to u > 7.05 and is non-dynamical. This distinction is based on the analysis of the behavior 
of the apparent horizon, which lies along the — 0 line. It changes its position along u 
in the dynamical region and becomes null as u ^ oc. Its location there, along u — 0.84, 
determines the event horizon location in the spacetime [119]. The Cauchy horizon is situated 
at the future null infinity, i.e., at u = oc. It is invisible within the covered spacetime domain, 
but its existence is confirmed by the fact that the r = const, lines tend towards constant u 
with increasing advanced time. The line r = 0 is non-singular up to u = 3.04 and becomes 
a spacelike singularity for bigger values of advanced time. 

One of physical quantities, which plays a significant role during the interpretation 
of the results is the quasi-local Hawking mass [120]. Its value calculated for the spherically 
symmetric spacetime with two coupled gauge fields and Pjj^ is provided by 


m (?i, v) 


+ + ‘^Q^ + T^ + aDMTQ \ 

2 V a2 r2 ) ' 


(3.4) 


It describes the mass contained within a sphere of a radius r{u, v) on a spacelike hypersurface 
containing the point (i^, v). The mass of a particular object is a value of the above expression 
at the event horizon in the non-dynamical region of the emerging spacetime, i.e., ior v — Vf. 


4 Electrically charged scalar field collapse with dark energy 

An interpretation of the dark energy impact on the examined collapse was made on the ba¬ 
sis of two distinct sets of solutions describing the outcomes of the process. The first one 
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corresponds to the varying amplitude of the phantom field pk and the constant amplitude 
of the electrically charged scalar field ps = 0.6. The other set was obtained during simula¬ 
tions with the constant phantom scalar field amplitude pk = 0.1 and the varying charged 
scalar field amplitude Such an approach allows to investigate the influence of the gravita¬ 
tional self-interaction strength of a particular field on the studied process. The evolutions 
conducted for the selected values of constant amplitudes, when the accompanying field 
is absent, lead to dynamical spacetimes of Reissner-Nordstrom and Schwarzschild types, 
respectively [16]. 

4.1 Dynamical emergence of wormhole-like structures and naked singularities 

Dynamical spacetimes resulting from the collapse proceeding for the constant electrically 
charged scalar field amplitude and varying pk are presented in figure 3. The emerging 
spacetimes are singular for small values of pk and become non-singular for its bigger values. 
This results from the fact that the phantom field contributes negatively to the overall energy 
of the system. For small values of the amplitude pk^ up to 0.22, the spacetimes are singular 
with both r^u = 0 and r^y = 0 lines of zero expansion visible in the domain of integration. 
For values of the altering amplitude not exceeding 0.185 there are two branches of the r^y = 0 
line. The outer one surrounds the singular part of r = 0 (for explanation, see section 3) 
and it becomes null as 'u ^ oc. The second branch of the line r^y = 0 lies inside the first 
one and its beginning coincides with the end of the singular spacelike r = 0 line. The line 
of r^u = 0 is situated beyond the second branch of r^y = 0 and begins at the same point. 
For small values of the parameter pk these two lines practically coincide. 

The structure described above, with two branches of r^y = 0 lines and the line r^u = 0 
beyond them, is a dynamical wormhole. The strong flare-out condition 0i^i > 0 with 
i — u^v [121] is fulfilled along the r^u — 0 line and the inner branch of r^y = 0. This 
fact proves that they are wormhole throats. Each of the throats is related to one travel 
direction. The outgoing throat is situated along r^u — 0 and the ingoing one along r^y — 0. 
The strong flare-out condition is not satisfied along the outer branch of the line = 0, 
because the relation 9y^y < 0 holds there. This confirms that the branch is not a throat, 
but an apparent horizon in the spacetime, whose null course in the region v ^ oo indicates 
the location of an event horizon. 

For values of the amplitude pk larger than 0.185, but not exceeding 0.22, the two 
branches of the line r^y = 0 form a loop, whose ends are joined to the ending points 
of the singular part of the r = 0 line. The line — 0 remains outside the loop. Such 
a structure corresponds to a naked singularity in the spacetime [9]. For values of the scalar 
field amplitude pk larger than 0.22, the spacetime is non-singular. 

The second set of solutions involves spacetimes emerging from a collapse when the am¬ 
plitude of the phantom scalar field is constant and the amplitude of the electrically charged 
scalar field changes. In figure 4 the structures of spacetimes which stem from the collapse 
are depicted. As in the previous set of solutions, the wormhole structures are observed. 
But unlike the preceding case, the spacetime is non-singular for small values of the ampli¬ 
tude ps^ not exceeding 0.35, and contains a wormhole for its bigger values. The wormhole 
structures are similar to these presented above, that is there are two branches of the r^y = 0 
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Figure 3. (color online) Penrose diagrams of spacetimes emerging from the SF-DE evolution 
with the varying phantom scalar field amplitude. The family parameter for the electrically charged 
scalar field Ps = 0-6, while for the phantom scalar field pk is equal to (a) 0.01, (b) 0.15, (c) 0.2 
and (d) 0.25. 


line, the outer and the inner one, and the line r^u — 0 located within them. The outer 
branch of the line r^y — 0 settles along a constant value of ix-coordinate for 'u ^ oo indicat¬ 
ing the location of an event horizon. Contrary to the previous case, no naked singularities 
were observed in the formed dynamical spacetimes. 

4.2 Vicinity of wormhole throats 

As was explained in the previous section, due to the behavior of the expansion (3.3) and its 
adequate derivatives, the ingoing and outgoing wormhole throats are located along the inner 
branch of r^y — 0 and along the — 0 line, respectively. As was anticipated and confirmed 
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Figure 4. (color online) Penrose diagrams 
of spacetimes emerging from the SF-DE evo¬ 
lution with the varying electrically charged 
scalar field amplitude. The family parameter 
of the phantom scalar field pk = 0.1, while Ps 
is equal to (a) 0.3, (b) 0.375 and (c) 0.7. 


for eternal wormholes, the violation of the null energy condition (NEC) 

> 0, (4.1) 

where denotes a null vector, is crucial for their existence [121, 122]. 

In spherical symmetry, taking into account the line element (2.12), the stress-energy 
tensor can be written in the general form 

Tj^iydx^dx^ = Tuudu^ + Tyydv^ + 2Tuvdudv + Tkk'^‘^dQ^^ (4.2) 

where Tyy, Tyy uud Tj^k its components, while = ©,$. The matter described 
by the above stress-energy tensor fulfills the null energy condition when diagonal compo¬ 
nents of (4.2) corresponding to the coordinates u and v are non-negative [123], i.e., 

Tuu > 0 and Tyy ^ 0. (4.3) 

The non-zero components of the stress-energy tensor (2.24) were calculated for the worm- 
hole spacetimes in order to examine the fulfillment of the NEC in the dynamical case. 
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Figure 5. (color online) The stress-energy tensor components (2.25)-(2.28) in the ('r’i4)-plane 
for evolutions carried out for electrically charged scalar field and phantom scalar field amplitudes 
set as (a) Ps = 0.6 and pk = 0.15, (b) Ps = 0.7 and pk = 0.1. 


Figure 5 presents the components (2.25)-(2.28) for the SF-DE collapse with the ampli¬ 
tudes of the electrically charged field and the phantom field equal to Ps — 0.6, pk = 0.15 
and Ps = 0.7, pk = 0.1 (the respective spacetime structures were presented in figures 3(b) 
and 4(c)). 

In the vicinity of wormhole throats the violation of the null energy condition is observed, 
because the requirement (4.3) is not satisfied. At the same time, this condition is fulfilled 
nearby the singular part of r = 0. The violation of NEC is thus a factor, which determines 
the existence of wormhole throats also in the dynamical case. Another area, in which 
the null energy condition is violated, is the vicinity of null lines v ^ 1.3 and u ^ 1.3 in both 
investigated wormhole spacetimes. This region is related to the evolution of the phantom 
scalar field, during which it neither falls into the central singularity nor evolves through 
the wormhole, but escapes to null infinity as u ^ oc. 

Figures 6 and 7 depict the modulus of the electrically charged complex scalar field 
and the phantom scalar field in the (u?i)-plane for the same evolutions as in figure 5. 
The modulus of the electrically charged scalar field increases with retarded time uniformly 
for all values of the u-coordinate. It is hard to distinguish any area, within which the field 
accumulates and the biggest values are one order of magnitude larger than the small¬ 
est ones. The value of the function related to the phantom field is considerably higher 
in the region nearby wormhole throats, where it is four orders of magnitude bigger than 
in the surrounding spacetime region. It can be interpreted as a phenomenon of stabiliza¬ 
tion of the wormhole structure by phantom matter [70]. The value of the phantom scalar 
field function is also significantly bigger along the hypersurfaces v ^ 1.3 and then ps 1.3 
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Figure 6. (color online) The (i;i4)-distribution of (a) the modulus of the electrically charged scalar 
field, |'0|, and (b) the phantom scalar field, 0, for the SF-DE evolution characterized by initial 
field amplitudes Ps = 0-b and pk = 0.15. 
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Figure 7. (color online) The (i;r4)-distribution of (a) \pj\ and (b) 0 for the SF-DE evolution 
characterized by initial field amplitudes ps = 0.7 and pk =0.1. 


indicating the course of the dispersive evolution of the field, which was also noticed during 
the above interpretation of the stress-energy tensor components. 

4.3 Black hole-wormhole duality 

A possibility of a two-way interconversion between black holes and wormholes was discussed 
in [124-128]. It was stated that a black hole may turn into a wormhole when irradiated 
by matter with negative energy. On the other hand, a wormhole stabilized by such exotic 
matter shall turn into a black hole when the matter source disappears. 

Dynamical wormhole spacetimes obtained in our computations contain spacelike singu¬ 
larities for large values of retarded time and the existence of event horizons is due to these 
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late-i/ singularities. A singularity exists during the dynamical stage of the collapse and bi¬ 
furcates into two wormhole throats when the spacetime tends towards its stationary phase 
at 'L’ ^ oc. Such a course of the investigated process suggests that a black hole, which 
emerges in the preliminary stage of the evolution, finally converts into a wormhole. The ex¬ 
istence of a singularity at large values of retarded time is strictly connected with the fact 
that the null energy condition is fulfilled in the area (figure 5), what is in turn a result 
of a considerably smaller value of the phantom field function in this region (figures 6(b) 
and 7(b)). The reason for the black hole-wormhole conversion is a substantial increase of its 
values in the particular area, which in consequence turns into a tunnel outlined by wormhole 
throats. 

4.4 Properties of wormholes 

The issue whether and in what manner do the properties of the obtained spacetimes differ 
when one of the field amplitudes is constant and the value of the other amplitude changes 
was also addressed. The ranges of varying amplitudes were chosen in such a way that 
in the case of the constant electrically charged scalar field amplitude the biggest pk is ap¬ 
proximately the last one not leading to the formation of a naked singularity. In the case 
of a constant phantom scalar field amplitude the smallest ps refers to the first dynamical 
wormhole spacetime which appears instead of a non-singular spacetime. 

In figure 8 the radii of event horizons and wormhole throats at the final value of the v- 
coordinate were plotted as functions of the altering field amplitudes. The chosen value 
of advanced time Vf corresponds to the well-established stationary phase of the evolution 
due to the fact that in all considered cases the apparent horizon r^y = 0 is a null hypersurface 
there. The common feature of both sets of solutions is that the radii of event horizons as well 
as wormhole throats increase when the varying amplitude increases. The radii of event 
horizons for the biggest values of pk are the only exception (figure 8(a)). In this region 
two branches of the line r^y = 0 approach each other, which for values of pk exceeding 
0.185 results in their connection and naked singularity formation. Such a tendency is not 
visible for the other set of solutions with the varying amplitude Ps (figure 8(b)), hence 
a naked singularity does not form. The changes of radii of event horizons are considerably 
more noticeable in the case of a varying amplitude of the electrically charged scalar field. 
On the contrary, the radii of wormhole throats vary more significantly when the phantom 
scalar field amplitude alters. The radii of both wormhole throats are identical for all values 
of the varying ps. For altering p^ they are equal up to some point and than they begin 
to differ more and more significantly. 

Figure 9 depicts the i^-locations of event horizons, wormhole throats at the final value 
of advanced time and the points, where the line r = 0 becomes singular (singularity 
origins) for both sets of solutions. When the phantom scalar field amplitude increases, 
the ^/-locations of all the mentioned features increase (figure 9(a)), while in the other case 
the opposite tendency is observed (figure 9(b)). Similarly to the radii, the ^/-locations 
of wormhole throats are exactly the same in the case of altering ps and they bifurcate 
at some point when pj^ increases. It is connected with the forthcoming naked singularity 
formation. The changes in ?i-locations are considerably bigger when the phantom scalar 
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Figure 8. The radii of event horizons and wormhole throats at the final value of the r^-coordinate 
as functions of (a) the phantom scalar field amplitude pk and (b) the electrically charged scalar 
field amplitude ps for the SF-DE collapse. The constant free family parameters for the particular 
evolutions were chosen to be equal to Ps =0.6 and pk = 0.1, respectively. 
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Figure 9. The locations of event horizons, wormhole throats at the final value of the i;-coordinate 
and the points, where the line r = 0 becomes singular, as functions of (a) pk and (b) Ps for the SF- 
DE collapse. The respective constant free family parameters are the same as in figure 8. 


field amplitude is constant and the amplitude of the electrically charged scalar field varies, 
in comparison to the other set of solutions. 

The masses of wormholes versus the phantom scalar field and electrically charged scalar 
field amplitudes are presented in figure 10. In both cases the mass increases when the am¬ 
plitude raises. 
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Figure 10. The wormhole masses, ttiehi as functions of (a) pk and (b) Ps for the SF-DE collapse. 
The respective constant free family parameters are the same as in figure 8. 

5 Influence of dark matter on electrically charged scalar fleld evolution 

During examining the role of dark matter in the studied collapse the values of parameters, 
which characterize the considered dark matter model, i.e., a^Mi and Xdm^ were var¬ 
ied within the permissible ranges outlined in section 2. It turned out that a^M and Xdm 
do not play a significant role during the process. The latter result is an extension of the pre¬ 
vious attempt to include the X^-type interactions in the gravitational collapse mentioned 
in [129], which was made for the evolution of a real scalar field without electric charge. 
The results of the conducted simulations will be presented for a^M — 10“^ and Xdm = 0.1. 
The course of the SF-DM evolution was also interpreted regarding the strength of a grav¬ 
itational self-interaction of the evolving fields by changing the values of their initial ampli¬ 
tudes ph and Ps. 

5.1 Spacetime structures 

The Penrose diagrams of spacetimes formed during the SF-DM collapse with a varying 
amplitude ph are presented in figures 11 and 12 for m? = 0.1 and the electrically charged 
scalar field amplitude Ps equal to 0.6 and 0.2, respectively. The structures of spacetimes, 
which stem from the process in the absence of dark energy are shown in figure 2. The former 
is the dynamical Reissner-Nordstrom spacetime and the latter is non-singular. 

In the case of ps — 0.6, all the resulting spacetimes are dynamical Reissner-Nordstrom 
spacetimes for all initial amplitudes of the complex field, which constitutes the dark mat¬ 
ter model. The central spacelike singularity at r = 0 is surrounded within each of them 
by a single apparent horizon situated along the = 0 line, which settles down along 
an event horizon as ^ oc. The Cauchy horizon is located at future null infinity v — oo. 
An increasing value of ph results in an earlier formation of a black hole in terms of retarded 
time. When ps equals 0.2, the obtained spacetime is non-singular for values of ph not ex¬ 
ceeding 0.135 or contains a dynamical Reissner-Nordstrom black hole for its larger values. 
This result is particularly interesting due to the role played by dark matter in this case. Al¬ 
though the self-interaction of an electrically charged scalar field is too weak to form a black 
hole (figure 2(a)), the presence of dark matter enables the emergence of a singular space- 
time. The appearance of the singularity is clearly related to the dark matter action during 
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Figure 11. (color online) Penrose diagrams of spacetimes formed during the collapse of the elec¬ 
trically charged scalar field of the amplitude Ps = 0.6 in the presence of dark matter with = 0.1 
and ph equal to (a) 10“^, (b) 0.1, (c) 0.15 and (d) 0.2. 
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Figure 12. (color online) Penrose diagrams 
of spacetimes formed during the SF-DM evo¬ 
lution with Ps = 0.2, m? = 0.1 and ph equal 
to (a) 0.1 and (b) 0.4. 


the collapse, but the existence of a Cauchy horizon is strictly connected with the charge 
associated with the scalar field coupled with the Maxwell field. 

The spacetimes formed during the discussed collapse for ps = 0.6, ph = 0.1 and several 
negative values of are shown in figure 13. For m? larger than —3.85 the dynamical 
Reissner-Nordstrom spacetimes form, while for its smaller values the Cauchy horizon is ab¬ 
sent in the spacetime and hence dynamical Schwarzschild black holes stem from the evolu¬ 
tion. In the latter case, the central r = 0 spacelike singularity is surrounded by an apparent 
horizon r^y = 0, which coincides with an event horizon at u ^ oo. The Cauchy horizon 
does not exist, because the r = const, lines do not settle along u = const, hypersurfaces 
as u ^ oo. For m? > 0 all the spacetimes are of the Reissner-Nordstrom type. The non-zero 
vev of the dark matter complex scalar field, i.e., the existence of dark matter in the form 
of a massive gauge boson (called Z' or dark photon) during the course of the studied process 
for < —3.85, results in the disappearance of a Cauchy horizon and favors the formation 
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Figure 13. (color online) Penrose diagrams of spacetimes formed during the SF-DM collapse 
with ps = 0.6, ph = 0.1 and equal to (a) —1, (b) —2, (c) —3 and (d) —4. 


of a simpler Schwarzschild-type spacetime. 

5.2 Fields behavior in spacetimes 

The components of the stress-energy tensor (2.24) calculated for the spacetimes formed 
during the SF-DM collapse characterized by the amplitudes Ps — 0.6, = 0.1 and m? 

equal to —1 and —4 are presented in figure 14 (the respective spacetime structures were 
presented in figures 13(a) and 13(d)). The ('Uix)-distributions of the moduli of the electrically 
charged scalar field and the dark matter complex scalar field for the same evolutions are 
shown in figures 15 and 16. 

Both Tuu and Tyy components of the stress-energy tensors are positive nearby the sin¬ 
gular r = 0 lines in the presented spacetimes. Hence, the null energy condition (4.1) 
is fulfilled in the vicinity of the central singularities. They do not bifurcate into worm- 
hole throats and black holes form during the SF-DM evolution. Similar conclusions stem 
from the analysis of the stress-energy tensor components for singular spacetimes presented 
in figures 11 and 12. 

An inspection of the moduli of the field functions reveals that the type of the forming 
spacetime depends on their values in the vicinity ofr^Oasr’^oo. A significant value 
of the modulus of an electrically charged scalar field in this region, three orders of magnitude 
bigger than in the remaining spacetime, accompanied by a small value of the other field mod¬ 
ulus, results in the formation of a dynamical Reissner-Nordstrom spacetime. It is connected 
with the repulsive character of the scalar field due to the electric charge associated with it. 
On the other hand, when the modulus of the dark matter complex scalar field is bigger 
at V ^ oo nearby r = 0 (two orders of magnitude in relation to the rest of the spacetime 
in the considered case), the importance of the other field is reduced and the dynamical 
Schwarzschild spacetime emerges. Such behavior of dark matter during the gravitational 
collapse is consistent with its action on a cosmological scale, where it tends to concentrate 
in confined areas due to the attractive character of its gravitational interaction. It is worth 
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Figure 14. (color online) The stress-energy tensor components (2.25)-(2.28) in the ('r’rt)-plane 
for SF-DM evolutions with ps = 0.6, ph =0.1 and equal to (a) —1 and (b) —4. 


(a) 


(b) 



Figure 15. (color online) The (i;i4)-distribu- 
tion of the moduli of (a) the electrically charged 
scalar field, It/^I, and (b) the dark matter com¬ 
plex scalar field, |x|, for the SF-DM evolution 
with Ps = 0.6, Ph = 0.1 and mS = —1. 


Figure 16. (color online) The ('r’rt)-distribu- 
tion of (a) and (b) |x| for the SF-DM evo¬ 
lution with Ps = 0.6, Ph =0.1 and mS = —4. 


noting, that in the regime of the selected dark matter model, the described behavior is re¬ 
lated to the change in potential (the value of m^), not in the gravitational self-interaction 
strength of the complex field (the amplitude ph)- 
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Figure 17. (color online) The ix-locations of the event horizons, ueh^ the masses, tueh^ 
and the radii, veh^ of black holes formed during the SF-DM collapse as functions of (a) ph 
with = 0.1 and (b) with = 0.1. The green line indicates the border between spacetimes 
of the Schwarzschild (m^ < —3.85) and Reissner-Nordstrom (m^ ^ —3.85) types. 


5.3 Characteristics of the formed black holes 

The dependence of the event horizon ?i-locations, the masses and the radii of black holes 
formed during the collapse of an electrically charged scalar field accompanied by dark matter 
on Ph and is depicted in figure 17. The maximal value of the initial amplitude of the dark 
matter scalar field was selected so that for all investigated cases the collapse begins outside 
of the event horizon. 

The black holes form earlier in terms of retarded time and their masses and radii 
increase as ph gets bigger (figure 17(a)). When m? increases the moment of the event 
horizon formation appears later in terms of u and the black hole masses and radii decrease 
with an inflection nearby m? — —3.85, which is the point of the change between the Reissner- 
Nordstrom and Schwarzschild spacetime structures (figure 17(b)). 

6 Collective effect of dark components on the considered collapse 

Similarly to the evolutions described in section 5, the collapse of an electrically charged 
scalar field with both dark energy and dark matter will be presented for the values of aE>M 
and Xdm equal to 10“^ and 0.1, respectively. The altering of the initial amplitudes 
Ph and pk is equivalent to changing the strength of the particular field gravitational self¬ 
interaction within the examined physical system. 

6.1 Spacetime structures 

The structures of spacetimes formed in the course of the SF-DE-DM collapse for ps — 0.6, 
Pk = 0.15, = 0.1 and the varying amplitude of the dark matter complex scalar field are 

presented in figure 18. The obtained spacetimes contain dynamical wormholes. The central 
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Figure 18. (color online) Penrose diagrams of spacetimes formed during the gravitational evolution 
of the electrically charged scalar field of the amplitude ps = 0.6 in the presence of the phantom 
scalar field with p;. = 0.15 and dark matter, whose model is described by the parameters = 0.1 
and ph equal to (a) 10“^, (b) 0.1, (c) 0.2 and (d) 0.25. 


spacelike r = 0 singularity bifurcates into two wormhole throats located along the lines 
r^y = 0 and = 0 as 'L’ ^ oc. The whole is surrounded by an apparent horizon = 0, 
which coincides with the event horizon in the region, where the spacetime settles down 
after the dynamical stage of the process. The value of ph does not influence the type 
of the emerging spacetime, but it causes the event horizon to form earlier in terms of retarded 
time. 

The Penrose diagrams of spacetimes which stem from the studied collapse with a vary¬ 
ing phantom field amplitude with p^^O.G, p/i = 0.15 and m? = 0.1 are shown in figure 19. 
The spacetimes contain either wormholes or naked singularities. The former emerge for pj^ 
smaller than 0.23 and the latter for its larger values. The i^-location of the event hori¬ 
zon increases with the value of the altering amplitude. The gravitational self-interaction 
strength of the phantom field decides about the type of a spacetime structure, which forms 
in the process. 

Figure 20 presents the results of the considered evolution for different negative val¬ 
ues of the parameter m^, when the initial field amplitudes are ps — 0.6, pk = 0.15 
and ph = 0.15. There exist dynamical wormholes in the emerging spacetimes for m? ^ —2.1 
and Schwarzschild black holes for its smaller values. The collapse with positive m? always 
results in a spacetime, which contains a dynamical wormhole. 

The above results indicate that when the field gravitational self-interaction strength 
is concerned, the role of dark energy coexists with the influence of dark matter. The dark 
energy component decides what type of an object forms. Dark matter does not interfere 
with the structure type, but regulates the amount of time which is needed for the in¬ 
trinsic object to form. When the impact of the dark matter potential is considered, 
it turns out that the value of m? influences the observed structure by controlling whether 
the wormhole throats remain open in the spacetime. Their absence results in the formation 
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Figure 19. (color online) Penrose diagrams of spacetimes formed during the collapse of the electri¬ 
cally charged scalar field of the amplitude Ps = 0.6 in the presence of dark matter with ph = 0.15, 
= 0.1 and the phantom scalar field with p^ equal to (a) 0.05, (b) 0.14 and (c) 0.25. 
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Figure 20. (color online) Penrose diagrams of spacetimes formed during the SF-DE-DM collapse 
with Ps = 0.6, pk = 0.15, ph = 0.15 and mn? equal to (a) —1, (b) —2, (c) —3 and (d) —5. 


of a Schwarzschild spacetime and the existence of a central singularity even after the dy¬ 
namical phase of the collapse is finished. It should be emphasized that the dark sector 
suppresses the natural tendency of an electrically charged scalar field to form a dynamical 
Reissner-Nordstrom spacetime in the course of the gravitational collapse. 

6.2 Fields behavior in spacetimes 

The (r’i/)-distributions of the stress-energy tensor components (2.25)-(2.28) for the dynami¬ 
cal spacetimes emerging from the SF-DE-DM evolutions characterized by the amplitudes 
Ps — 0.6, Pk = 0.15, Ph = 0.15 and m? equal to —1 and —3 are presented in figure 21 (the re¬ 
spective spacetime structures were presented in figures 20(a) and 20(c)). The spacetime 
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(a) 



Figure 21. (color online) The stress-energy tensor components (2.25)-(2.28) in the ('r’rt)-plane 
for SF-DE-DM evolutions with ps = 0-6, Pk = 0.15, ph = 0.15 and equal to (a) —1 and (b) —3. 


distributions of the moduli of the electrically charged scalar field, the dark matter complex 
scalar field and the phantom scalar field for the same cases are shown in figures 22 and 23. 

The Tuu and Tyy components of the stress-energy tensors are negative in the vicinity 
of wormhole throats when m? = —1. The null energy condition (4.1) is violated there. 
In the case of m? = — 3 the components are positive nearby r = 0 and hence a black 
hole, not a dynamical wormhole arises in the spacetime. The obtained result strengthens 
the conclusion presented in section 4.2 that the violation of NEC is needed for the formation 
of wormholes during dynamical evolutions. Clearly, this statement is valid even for compli¬ 
cated matter-geometry systems evolutions, in which several types of interacting matter are 
involved. 

Similarly to the SF-DM case, the obtained spacetime type depends on the values 
of the evolving field functions in the vicinity of r = 0. The formation of a dynamical worm- 
hole stems from a significant value of the modulus of the phantom scalar field in this region 
as u ^ oo (four orders of magnitude bigger than in the remaining spacetime) accompanied 
by a small value of the dark matter field modulus. When the modulus of the dark matter 
complex scalar field is bigger at u ^ oo nearby r = 0 (three orders of magnitude in rela¬ 
tion to the rest of the spacetime in the considered case), the importance of the remaining 
fields is reduced and the dynamical Schwarzschild spacetime emerges. The structure type 
and the fields behavior in a spacetime depend on the value of the parameter m^, not the ini¬ 
tial amplitude of either of the collapsing fields. 
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Figure 22. (color online) The ('r’i4)-distribution of (a) the moduli of the electrically charged scalar 
field, |'0|, (b) the moduli of the dark matter complex scalar field, |x|, and (c) the phantom scalar 
field, 0, for the SF-DE-DM evolution with ps = 0.6, pk = 0.15, ph = 0.15 and = — 1. 



Figure 23. (color online) The (i;i4)-distribution of (a) \'ip\, (b) |x| and (c) (j) for the SF-DE-DM 
evolution with ps = 0.6, pk = 0.15, ph = 0.15 and = —3. 


6.3 Properties of the emerging objects 

The influence of ph^ pk and on the event horizon ^/-locations, the masses and the radii 
of objects formed during the gravitational evolution of an electrically charged scalar field 
accompanied by dark matter and dark energy is shown in figure 24. The maximal value 
of the amplitude ph was determined by the event horizon location in the spacetime so as to be¬ 
gin the collapse outside it. The maximal considered value of the phantom scalar field ampli¬ 
tude was the biggest value not leading to the formation of a naked singularity in spacetime. 

As the dark matter complex scalar field amplitude increases, the wormholes appear ear¬ 
lier in their radii increase and masses decrease up to = 0.23 and increase for larger ph 
(figure 24(a)). In the case of an increasing pk^ the values describing all the proper¬ 
ties increase, apart from veh for values of the amplitude exceeding 0.21 (figure 24(b)). 
Just as in the case of the SF-DE system, such behavior indicates the upcoming forma- 
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Figure 24. (color online) The i^-locations 
of the event horizons, uehi tbe masses, rriEH^ 
and the radii, veh^ of the objects formed 
during the SF-DE-DM collapse as functions 
of (a) Ph with Pk = 0.15, = 0.1, (b) pk 

with Ph = 0.15, w? = 0.1 and (c) w? with 
Pk = 0.15, Ph = 0.15. The green line indicates 
the border between spacetimes containing black 
holes of the Schwarzschild type {m? < —2.1) 
and wormholes {w? ^ —2.1). 


tion of a naked singularity in the spacetime. When m? increases the event horizon forms 
later in terms of retarded time and the wormhole radii decrease with an inflection nearby 
m? — —2.1, which is the point of the change between the wormhole and Schwarzschild 
spacetime structures (figure 24(c)). The behavior of masses is more complicated, as apart 
from the inflection around m? = —2.1, there exists a minimal mass for m? — —3.9. 
The non-monotonic behavior of the described characteristics is related to the fact that 
the investigated system consists of multiple components and the observed phenomena are 
a compromise among the tendencies characteristic for particular evolving matter types. 

7 Conclusions 

The influence of dark matter and dark energy on the electrically charged scalar field collapse 
was investigated. Dark matter was described by a complex scalar field with a quartic 
potential, coupled to a (7(l)-gauge field. Dark energy was modeled by a scalar field coupled 
to gravity in a phantom manner. The exponential relation of the scalar field with dark 
energy originated from the low-energy limit of the string theory. Dark matter was coupled 
with the electrically charged scalar field through a kinetic coupling between the two gauge 
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fields present in the system. 

The existence of phantom matter in spacetime can result in the formation of naked sin¬ 
gularities or dynamical wormholes during the gravitational collapse of an electrically charged 
scalar field. The latter objects have two throats related to two travel directions and are 
surrounded by event horizons. The values of the phantom field are considerably bigger 
in the spacetime region situated nearby wormhole throats, in comparison to the surrounding 
area. This accumulation leads to the violation of the null energy condition, which is essen¬ 
tial for the dynamical wormhole formation and the stabilization of its throats. The masses 
of wormholes, the radii of event horizons and wormhole throats increase when the ampli¬ 
tudes of both phantom and electrically charged scalar fields increase. The event horizon, 
wormhole throats and singularity origin appear later in terms of retarded time as pk in¬ 
creases while ps is constant and at smaller values of the ^/-coordinate in the opposite case. 

In the obtained wormhole spacetimes spacelike singularities situated along r — 0 exist 
during the dynamical stage of the evolution. They turn into wormhole throats as the space- 
time tends towards the final stage of the collapse. This suggests that the considered space- 
time shortcuts ensue from the process, during which a black hole gains exotic matter of nega¬ 
tive energy. Thus the existence of duality between black holes and wormholes was confirmed 
in the course of fully non-linear computations. 

The spacetimes resulting from the examined collapse in the presence of dark matter are 
either non-singular or contain black holes, depending on the gravitational self-interaction 
strengths of both evolving fields, that is their initial amplitudes Ps and p^. The type 
of an emerging object depends on the value of the square of mass parameter from the dark 
matter model. In general, the non-zero vev of the dark matter complex scalar field favors 
the formation of dynamical Schwarzschild spacetimes for w? < —3.85, while for the vanish¬ 
ing vev the dynamical Reissner-Nordstrom black holes appear. 

The null energy condition is fulfilled in the vicinity of singularities in the emerging 
spacetimes, hence the stabilization of wormhole throats is impossible and such objects 
do not form. The values of the moduli of the electrically charged scalar field and the dark 
matter complex scalar field are considerably bigger nearby r = 0, where ^ oo, than 
in the remaining regions in spacetimes of Reissner-Nordstrom and Schwarzschild types, 
respectively. This indicates a repulsive character of the former and an attractive character 
of the latter field. 

The moment of a black hole formation, i.e., an appearance of an event horizon in terms 
of retarded time, as well as the black hole masses and radii are controlled by the values of ph 
and m?. The event horizon //-locations increase and the remaining characteristics decrease 
as the dark matter complex scalar field initial amplitude raises. The opposite tendency 
is observed in the case of an increasing parameter m?. 

The dark sector suppresses the tendency of a self-interacting electrically charged scalar 
field to form a dynamical Reissner-Nordstrom black hole during the gravitational evolution. 
In the case of its overall influence on the investigated process, the strength of the dark en¬ 
ergy gravitational self-interaction prevails over the dark matter gravitational self-interaction 
strength. The dark energy component decides on the type of an object, while dark matter 
regulates the amount of time which is needed for it to form. The parameter m? of the dark 


-32- 



matter potential influences the observed structures by controlling whether the wormhole 
throats appear in the spacetime. Their absence results in the formation of a Schwarzschild- 
type black hole, for the non-zero vev of the dark matter complex scalar held and m? < —2.1. 

In the spacetimes resulting from the SF-DE-DM collapse the null energy condition 
is violated in the vicinity of wormhole throats and fulfllled nearby the singularities, just 
as in the previous cases. The formation of a wormhole requires values of a phantom scalar 
held considerably bigger nearby the throats in comparison to the rest of a spacetime. Large 
values of the modulus of the dark matter complex scalar held result in the emergence 
of a Schwarzschild-type object. 

The behavior of ?i-locations, masses and radii as functions of ph^ pk and m? are com¬ 
plicated due to the fact that the investigated system consists of multiple components, each 
of which displays a different tendency during the evolution. The electrically charged scalar 
held favors the Reissner-Nordstrom black hole formation, dark energy supports the emer¬ 
gence of wormholes or naked singularities and dark matter tends to form the Schwarzschild- 
type objects. The observed phenomena are a compromise among these tendencies charac¬ 
teristic for the involved matter types. 

The changes introduced by the dark sector to the electrically charged scalar held grav¬ 
itational collapse are presented in table 1. It summarizes the obtained results in the con¬ 
text of the role of particular parameters in the studied evolution for the initial amplitude 
of an electrically charged scalar fleld, which solely leads to the formation of a black hole. 
The obtained outcomes confirm the previous conclusions regarding the role of dark energy 
during the gravitational collapse. It does not prevent the formation of singular spacetimes 
during the dynamical evolution when it is accompanied by other types of matter, in partic¬ 
ular dark matter [18, 20]. The current studies also widened previous findings in this regard. 
They allowed us to describe the structures of objects existing in these singular spacetimes, 
which turned out to be wormholes, naked singularities or black holes of Schwarzschild 
or Reissner-Nordstrom types. 

The Penrose diagrams of spacetimes, on which the spacetime structures were presented 
in the current paper, are related to the Carter-Penrose diagrams via conformal transfor¬ 
mation, which preserves the causal structure of a spacetime. The Carter-Penrose diagrams 
of both static and dynamical Schwarzschild and Reissner-Nordstrom spacetimes can be 
found, e.g., in [6, 8, 1 14]. In figure 25 we present such diagrams for the dynamical wormhole 
and naked singularity spacetimes, which stem from the performed analyses of the dynamical 
formation of these objects. 

Apart from investigations of spacetime structures formed during the collapse in the pres¬ 
ence of the dark sector, other tantalizing questions are the mechanism of black hole forma¬ 
tion and the behavior of matter accreting onto a nascent black hole taking the cosmological 
evolution into account. It was claimed that primordial black holes can accrete the surround¬ 
ing dark energy, phantom energy or ghost condensate very effectively [130, I3I]. However, 
the conducted analyses neglected cosmic expansion. 

Until now, there have been several attempts to And black hole solutions embedded 
in the expanding Universe, starting from the Einstein-Straus [132] and McVittie [133] solu¬ 
tions to the multi black hole Kastor-Traschen cosmological solution [134], which is a time- 


- 33 - 


Evolution 

Pk /' 

Ph Z’ 

m? 

SF-DE 

WH ^ NS ^ ns [16] 
Ueh E 

Teh Z' uieh Z' 

- 

- 

SF-DM 

- 

RN 

Ueh \ 

rEH Z' uieh Z' 

RN ^ S 

Ueh Z' 

rEH \ ruEH \ 

SF-DE-DM 

WH NS 

Ueh Z' 

rEH E\ iueh Z' 

WH 

Ueh \ 

rEH Z’ ruEH \Z' 

WH S 

Ueh Z' 

rEH \ iueh \Z' 


Table 1. Spacetime structures formed during the gravitational collapse of an electrically charged 
scalar field with Ps = 0.6 accompanied by dark energy and dark matter and the selected charac¬ 
teristics of the emerging objects. Symbols: ns - non-singular spacetime, S - Schwarzschild-type 
spacetime, RN - Reissner-Nordstrom-type spacetime, WH - wormhole, NS - naked singularity, 
ueh - r^-location of an event horizon, teh and rriEH - radius and mass of the intrinsic object, 
respectively. Arrows: /^/\ - an increase/decrease of a value of the particular characteristic or am¬ 
plitude, ^ - the change of the spacetime structure while the particular parameter varies. 



Figure 25. (color online) The Carter-Penrose diagrams of dynamical spacetimes, which contain 
(a) a wormhole and (b) a naked singularity. The red and blue solid lines are the = 0 and = 0 
horizons, respectively. The red and blue dashed line symbolizes wormhole throats. The meaning 
of the symbols is the same as in figure 1. The causal structure of the spacetime regions marked gray 
on the diagrams is impossible to interpret without any doubts due to the limitations of numerical 
computations in these areas. 

dependent generalization of the Majumdar-Papapetrou one. An exact analytical black hole 
solution in the expanding Universe was derived within the Einstein-scalar-Maxwell system 
with two U(l)-gauge fields and an exponential potential of the scalar field [135]. The dila- 
ton black hole on a thick brane and a cosmological brane black hole solutions were studied 
in [136, 137], while the time-dependent solution from compactification of intersecting branes 
in higher-dimensional unified theories was found in [138]. As was already stated in the In- 
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troduction, the dynamical solutions may behave quite differently in comparison to their 
stationary cousins. Some realistic attempts to analytically describe the accretion of a mat¬ 
ter field onto an evolving black hole were conducted in [139], where a time-dependent scalar 
field and a time-dependent black hole were considered. In order to correctly account for 
the scalar field accretion, the expansion in terms of a slow roll parameter was implemented. 
The only initial assumptions about the geometry were its SO (3) symmetry and the time 
and radial dependence of the line element. Hence, the presented method did not require 
a direct assumption about the form of the underlying black hole metric. 

The issues introduced above to be studied in the fully non-linear dynamical context are 
beyond of the scope of the presented research, as they require different methodology than 
the one employed in the paper. We hope to return to the announced problems in relation 
to the cosmological evolution of the Universe in the future researches. 

A Numerical computations 

The numerical simulations of the investigated process were conducted with the use of a mod¬ 
ified version of the code, the specifics of which were presented in [14]. The scheme was broad¬ 
ened and adapted for the presently considered problem by including the phantom constant ^ 
and the additional sector (2.52)-(2.55), which refers to dark matter. Due to the coupling 
between the and fields, the module covering the Maxwell field evolution (2.51) was 
also modified. The module governing the geometry dynamics (2.44)-(2.47) was adjusted 
to the presently examined physical system. 

The inclusion of the additional sector in the code required posing adequate initial 
and boundary conditions for functions zi, Z 2 , /ii, /i 2 , '^ 2 , 7 , and a^ 2 - The first two 

evolved along u according to equations and The dynamics of /ii, /i 2 , wi^ W 2 ^ 

T and 7 along v was governed by equations P7, P 8 , respectively. 

The values of the remaining two quantities were calculated on the basis of their definitions 
included in (2.35). 

Initial conditions were posed on a null hypersurface denoted as = 0. The profiles 
of functions hi and /i 2 were assigned according to (3.2). The values of and Z 2 were calcu¬ 
lated analytically using the relations P7 and P8. The behavior of the rest of the quantities, 
i.e., wi, W 2 ^ T and 7 , at the initial spacetime slice was resolved with the use of the three- 
point Simpson method and the Newton method at the first point. 

The line u — v was the line of boundary conditions. Because of the fact that the col¬ 
lapsing physical system consists of a matter shell, the values of T and 7 were equal to zero 
at the boundary, which refers to the non-singular r = 0 line. The quantities tzi and a ^2 also 
vanished due to the equations and . The flattening of the functions hi and /i 2 

was imposed on the boundary, i.e., hiy and /i 2 ,r were equal to zero there. The values 
of the quantities wi and W 2 were the same as zi and 2 : 2 , respectively. 

The accuracy of the numerical code was checked indirectly, due to the fact that no an¬ 
alytical solutions exist for the investigated process. The tests were performed for the evo¬ 
lutions initiated with parameters ps = 0.6, ph — 0.15, = 0.1, pk = 0.14 and ps = 0.6, 

Ph = 0.15, m? = —3, Pk = 0.15. The respective spacetime structures are presented in fig- 
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Figure 26. (color online) The convergence of field functions. The phantom field, 0, and the moduli 
of complex scalar fields, \ 2 p\ and |x|, were plotted versus v for evolutions conducted with integration 
steps, which were multiples of h = 10“^, along hypersurfaces of constant u equal to (a) 1.3504 
for evolution 1 and (b) 1.2504 for evolution 2. 


ures 19(b) and 20(c). For the purpose of testing the numerical code, these two cases will 
be referred to as evolutions 1 and 2, respectively. The analyses were done for spacetimes 
obtained within the SF-DE-DM physical system, which was the most comprehensive 
among those discussed in the paper. Hence, the performed code accuracy checks covered 
the operation of all the constituent modules of the code. 

The first test was based on checking the convergence of the obtained results. In order 
to monitor the convergence, the computations for evolutions 1 and 2 were conducted on four 
grids with integration steps being multiples of 10“^. A step of a particular grid was twice 
the size of a denser one. The convergence was examined onu = const, hypersurfaces chosen 
for each case arbitrarily. The selected hypersurfaces were situated close to the emerging 
horizon, but in the region where the adaptive grid was yet inactive, which was necessary 
for performing a proper comparison of the outcomes. 

The field functions along the selected hypersurfaces of constant u are shown in fig¬ 
ure 26. The areas in which the differences among functions obtained on various grids were 
most significant were magnified. The maximum observed discrepancy between the finest 
and coarsest grids was equal to 1.5%. Figure 27 presents the linear convergence of the nu¬ 
merical code. The maximal divergence between the field profiles obtained on two grids 
with a quotient of integration steps equal to 2 and their respective doubles was 5.7%. 
As expected, the errors became smaller linearly as the grid density increased. The over¬ 
all convergence analysis revealed that both the algorithm and the numerical code were 
adequate for solving the system of equations (2.36)-(2.55). 

The second test of the numerical code was based on checking the mass and charge con- 
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Figure 27. (color online) The linear convergence of the code. The differences between the phantom 
field functions, A(/), and the moduli of complex scalar fields, A|'0| and A|x|, calculated on grids with 
different integration steps (multiples oi h = 10“^) and their doubles were obtained along the same 
hypersurfaces of constant u as in figure 26 for (a) evolution 1 and (b) evolution 2. 


servation in the spacetime. The Hawking mass (3.4) and the charges related to both complex 
scalar fields coupled to the C/(l)-gauge fields (2.14) and (2.17) as functions of retarded time 
along the line 'L’ = 7.5, which was a maximal value of advanced time achieved numerically, 
are presented in figure 28 for the investigated evolutions. Since during the process the mat¬ 
ter was scattered by the gravitational potential barrier when the collapsing shell approached 
its gravitational radius, the physical quantities were not conserved during the whole evo¬ 
lution. The effect of the outgoing flux was negligible except for the vicinity of the event 
horizon. The deviation from the constancy increased with advanced time, as the horizon 
was approached. The maximal percentage deviations from the particular quantity conser¬ 
vation up to = 1.25 were equal to 3.3%, 8.1% and 1.1% for the mass, the electric charge 
and the charge associated with the dark matter sector, respectively. The analysis of mass 
and charge conservation in spacetime led to the conclusion that the behavior of matter 
investigated numerically was correct within the domain of integration. 
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Figure 28. (color online) Hawking masses, m/, and charges related to the [/(l)-gauge fields, Qf 
and Tf, calculated along i;/ = 7.5 as functions of retarded time, for both tested evolutions. 
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